arXiv:1503.06756v2 [cond-mat.mtrl-sci] 15 Sep 2015 


Out of equilibrium anomalous elastic response of a water nano-meniscus 


Simon Carpentier/’^ Mario S. Rodrigues,^ Miguel V.Vitorino,^ 

Luca Costa,^ Elisabeth Charlaix,^’^ and Joel Chevrier^’^ 

^ Univ. Grenoble Alpes, F-38000 Grenoble, France 
^ GNRS, Inst NEEL, F-38042 Grenoble, France 

^ Uniservity of Lisboa, Falculty of Sciences, BioISI-Biosystems & Integrative Sciences Institute, Gampo Grande, Lisboa, Portugal 
‘^ESRF, The European Synchrotron, 71 Rue des Martyrs, 38000 Grenoble, Erance 
^GNRS, LIPhy, Grensoble, F-38402, France 

We report the observation of a transition in the dynamical properties of water nano-menicus 
which dramatically change when probed at different time scales. Using a AFM mode that we name 
Force Feedback Microscopy, we observe this change in the simultaneous measurements, at different 
frequencies, of the stiffness G’(N/m), the dissipative coefficient G”(kg/sec) together with the static 
force. At low frequency we observe a negative stiffness as expected for capillary forces. As the 
measuring time approaches the microsecond, the dynamic response exhibits a transition toward 
a very large positive stiffness. When evaporation and condensation gradually lose efficiency, the 
contact line progressively becomes immobile. This transition is essentially controlled by variations 
of Laplace pressure. 


Visco-elastic properties of water nanobridges[T] at very 
different time scales, have never been investigated despite 
ubiquitous presence of capillarity. Associated forces are 
among the most intense at nanoscales with important 
consequences in soils and granular media. Interest in 
dynamical properties is immediately raised if one con¬ 
siders interacting surfaces with roughness scales down to 
nanometer. Even at moderate speeds, such as v=lm/s, 
characteristic times of surface interaction down to mi¬ 
crosecond appear in these conditions. Our measurements 
approaching these time scales, further strengthen the rel¬ 
evance of the dynamical properties to describe how real 
surfaces interact and are certainly of crucial importance 
in numerous AFM experiments |2] . We here report mea¬ 
surements of dynamical properties of a water nanobridge 
for a continuous range of the surface gap and a frequency 
bandwidth up to 0.1 MHz. We identify two regimes: one 
is the thermodynamical equilibrium; the second is out 
of equilibrium. Evaporation and condensation of water 
molecules between the liquid and the gas phase ensures 
that the nano-meniscus curvature is the one at thermo¬ 
dynamical equilibrium(2H=-l/rk) where 2H is the wa¬ 
ter bridge curvature and rk is the Kelvin radius. At 
time short enough, molecule exchanges between the liq¬ 
uid and the gas phase are no longer efficient and the wa¬ 
ter nanobridge is led to acquire a constant volume. The 
liquid bridge relaxation time is the time needed for the 
bridge to adapt its shape as required by thermodynam¬ 
ical equilibrium, when its length h is abruptly changed 
by 8h. This is controlled by molecular transport through 
diffusion mechanisms in gas phase. This relaxation time 
T can be estimated as in SEA context, see Ref.jS]- 

T = 2-fpr‘^ln{R/p)/Psat D (1) 

D is the diffusion coefficient of water molecules in air 
(D= 0.282 cm^/sec from Ref.|l]). p is the bridge az¬ 
imuthal radius, r, its meridional radius. Pgat is the satu¬ 
rated vapor pressure, R is the typical scale of the overall 


system and y the surface tension of the vapor liquid in¬ 
terface. In our experimental conditions, with a Kelvin ra¬ 
dius rk of about 12nm, the characteristic time is found to 
be around t= 10“® second. As w/27i increases from 300Hz 
to 114kHz, in our experiment, cox evolves from 2.10'^ 
to 0.7. At WT 2.10“^, the liquid nanobridge should be 
observed at thermodynamical equilibrium. At cot= 0.7, 
w/2'n:=114kHz, this time x is too short for the thermo¬ 
dynamical equilibrium to settle during oscillation. We 
shall therefore consider an extreme regime with constant 
volume when dealing with analysis of experiments per¬ 
formed at these high frequencies. 

Using the AFM mode that we name Force Feedback 
Microscopy (FFM), we have measured simultaneously 
the static force, the interaction stiffness G’ and the as¬ 
sociated dissipation G” at different frequencies. As de¬ 
scribed in Fig[^(a), the FFM used to perform this ex¬ 
periment is based on a homemade AFM setup with op¬ 
tical fiber |^. These three simultaneous measurements 
are done using a single nanotip-microlever system, with 
excitation frequency varied from 300Hz up to 114kHz. 
During FFM experiments, a feedback force is applied in 
real time to the tip so that it cancels the force due to 
the surface (this suppresses the “jump to contact” ob¬ 
served in classical AFM force-approach curve). The to¬ 
tal force applied to the tip is then equal to zero and the 
tip remains immobile. To apply this force in real time, 
a piezoelement simply changes the DG position of the 
clamped end part of the microlever. This measured dis¬ 
placement multiplied by the lever stiffness k results in 
the static tip/surface force measurement. Measured on 
this basis HE], the static force is reported in FigJ^(b). 
Both the tip radius and the Kelvin radius, rk are then 
estimated. The minimum force is roughly -diiyR, which 
leads to R close to 14nm, see more details in Ref.[S]. The 
slope of the force versus distance as the AFM tip is re¬ 
tracted, gives rk=12nm nunj. As shown in Fig0(a), a 
nanometeric oscillation is imposed to the tip. Measure of 
the tip amplitude and the associated phase shift leads to 
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Figure 1. (a) Description of the FFM measurement strategy. 
The DC tip position is a constant. This is due to real time 
application of a feedback force determined through a control 
loop after measurement of the tip position. Piezodither en¬ 
forces an AC tip oscillation with 1 nm of free oscillation am¬ 
plitude. Aosc and $ are the measured quantities, (b), Static 
force measured between a brand new silicon AFM tip and a 
plane silicon sample during approach (blue) and retract (red) 
by FFM in presence of a capillary bridge. The linear behavior 
is htted using the approximate solution given by J.Crassous 
et al. in Ref. [8] with R=14nm and rk=12nm. 


the determination of the stiffness interaction and the as¬ 
sociated dissipation following these linear equations (see 

details in Ref. my 

G' = Fr[ncos($) — cos($oo)] (2) 


G” = — nsm(ti)) -I- sm($oo)] (3) 

UJ 

$ is the phase shift between excitation at the lever 
clamped end and the tip oscillation, n is the normalized 
amplitude of tip vibration i.e. the measured amplitude at 



Figure 2. From blue to orange curves (respectively 300 Hz, 1 
kHz, 44 kHz, 94 kHz and 114 kHz) show interaction stiffness, 
G’ versus tip surface distance as deduced from experimental 
measurements using EqlD Red curve is the negative of the 
numerical derivative of the measured static force as the tip is 
pulled away from the surface. Inset: G”(kg/sec) result from 
Eq|^ It characterizes the mechanical energy dissipation in¬ 
crease during tip oscillation due to the capillary bridge. G” at 
high frequencies is distinctively decreased as compared to G” 
at low frequencies. Dissipation when in the thermodynamical 
equilibrium is clearly much higher. 


distance z divided by the measured amplitude when tip 
is far away from the surface. At large distance n=l. In 
the experimental environment used (normal conditions 
in air), values of Fr and $oo are directly obtained from 
the system transfer function when the tip is far from the 
surface. 


F, = [(fc-mw2)2 + cV]i/2 (4) 


$00 = arctan[{c/m)w / (5) 

In the used frequency range with first resonance at 
74,330kHz, a single mode description of the lever dynam¬ 
ics is sufficient. The stiffness of the lever k is obtained 
from measurement of the Brownian motion. The damp¬ 
ing coefficient c and coq (and therefore m) are obtained 
from lever transfer function measurement when the tip is 
far from the surface. F,- and 4>oo determined using this 
method are inserted in Eqj^ and to obtain G’ and G” 
from n and <i>. In Figj^ k=2 N/m, c=2.49 10“® kg/sec 
and mo/27i=74,330kHz. 

Figj^ shows that the visco-elastic coefficients exhibit a 
strong variation as the excitation frequency is varied over 
close to 3 orders of magnitude: G’ varies from -IN/m up 
to lON/m (therefore from negative to positive). The pos¬ 
itive stiffness measured at high frequency, G’=+10N/m 
is a very high interaction stiffness in AFM field. It is sur¬ 
prisingly observed associated to water and to an intense 
attractive static force close to lOnN. From this apparent 
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positive stiffness, although only surface effects due to the 
nano-meniscus are at work, an effective Young modulus 
E can be formally deduced: E = G’z/n rb- It is in the 
range 0.1-1 GigaPa and increases as the tip moves closer 
to the surface. 

At o/2ii=300Hz, the measured interaction stiffness G’ 
corresponds to the numerical derivative of the static force 
versus distance (i.e. the derivative of the red curve in 
Figj^b)). Here wt<< 1 and the water nanobridge ap¬ 
pears in thermodynamical equilibrium. However, the 
gradual increase of the excitation frequency leads to an 
important increase of the interaction stiffness, which fi¬ 
nally becomes positive at all investigated distances m- 
The high frequency regime with positive and high stiff¬ 
ness is therefore an out of equilibrium behavior of the 
water droplet. As determined by use of Eq. an esti¬ 
mate of the bridge relaxation time x is 10“®sec: the water 
molecule exchange between liquid and gas diminishes as 
the frequency increases toward IMHz. This ultimately 
leads to a regime with a water nanobridge that has a 
constant volume. With this hypothesis of constant vol¬ 
ume at high frequency, we analyze the observed positive 
stiffness G’ due to the water nano-menicus properties. 
This in fact follows analysis of Lambert et al in Ref. m 
produced for far larger water bridge (millimeter range). 
Equation]^ shows the expression of the total force due to 
the water nano-menicus on the tip: 


F = 27rrb7 sin(0) — 2HTTr'^j (6) 

The force decomposes in two terms. The first term 
of Eq|^ is due to the tension associated to the contact 
line whereas the second is due to the Laplace pressure, 
y is the water surface tension, 2iirb the length of the cir¬ 
cular contact line (or triple line), (2H=l/p-l/r) and h, 
the angle between the water bridge and the surfaces. At 
thermodynamic equilibrium (u)t<<1) with efficient evap¬ 
oration condensation processes, whatever the tip surface 
distance and the oscillation amplitude, H is a negative 
constant. The pressure inside the bridge then does not 
change, and the contact line moves to accommodate this 
constrain of constant curvature as the tip surface distance 
is varied. During shortening of the bridge (8h<0), the ra¬ 
dius rb increases as a result of the bridge spreading: the 
resulting stiffness 8F/5h is negative as observed experi¬ 
mentally. At U)T>>1, we consider that the nanobridge 
volume remains constant. If the contact line does not 
move, the radius rb is constant and the curvature in¬ 
creases (H<0) as described in Figj^a). The second term 
in the force then provides a contribution to stiffness that 
is: 


SF 

Jh 


—27rr^7 


Jh 


( 7 ) 


This 8F/8h is then a positive contribution, which ex¬ 
plains the origin of the positive stiffness experimentally 
observed. 




Figure 3. (a), evolution of the capillary bridge shape as 

the gap decreases at constant volume and at locked con¬ 
tact line, (b), results from numerical calculation, reproduce 
the key experimental features obtained using Force Feedback 
Microscopy at different frequencies. At low frequency (blue 
curve), the signature of the thermodynamical equilibrium is 
again obtained with a negative stiffness at all distances. At 
high frequency oscillation (green curve), the water nanobridge 
is at constant volume and the contact line is immobile. A 
positive stiffness is obtained like in experimental results. It is 
however significantly above all measured stiffnesses for which, 
at cot= 0.7, contact line mobility is reduced but certainly far 
from zero. The orange curve represents a regime at constant 
volume, but with a contact line free to move. It leads to 
a negative stiffness. This regime, which implies a curvature 
only marginally varying, cannot explain the experimental ob¬ 
servations at high frequency. 


A numerical calculation of the stiffness G’ in the 
three extreme regimes (thermodynamical equilibrium, 
constant volume with locked contact line or free to move 
contact line) is reported in Figj^b). 

Finally as shown in inset of Fig|^ G”(kg/sec), the 
dissipative part of the linear response, is found to 
decrease as the excitation frequency increases and as 
G’ severely increases. This leads us to conclude that 
this change of G” is also determined by the contact 
line behavior. A mobile contact line on the AFM tip 
then appears to be associated to an important dissi¬ 
pation. As mobility of the contact line decreases, this 
channel of energy dissipation closes and G” diminishes 
accordingly. Although only surface effects are here 
at work, G” can be formally analyzed as an apparent 
volume viscosity: G”ocr)R, where R=14nm is the tip 
size. In that case, with a typical G”=10“® kg/sec here 
measured and taken equal to rjR, r) is in the range 
lO^-lO^Pa.sec. This is orders of magnitude above /)= 
lO'^Pa.sec, the bulk value of water viscosity (Ref. [HI [H]). 

In conclusion, for the first time, two different visco¬ 
elastic regimes in dynamical properties of a water 
nanobridge are measured. The major result is that, de¬ 
spite the fact that the static interaction between sur¬ 
faces associated to capillary bridges is always attrac¬ 
tive, an apparent and highly positive stiffness appears 
at the nanoscale when the two surfaces interact with a 
short characteristic time, whereas, the friction measured 
in same conditions, is found strongly reduced. A de¬ 
scription can be proposed to simultaneously explain these 
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results. Contrary to the situation at thermodynamical 
equilibrium with condensation and evaporation, in the 
out of equilibrium regime, any displacement of the line 
can only be ensured by a flow in the nanodroplet. The 
sticky boundary condition leads to a diverging shear ve¬ 
locity as the water film thickness goes to zero close to the 
triple line. In this situation, moving the line becomes ex¬ 
tremely difficult. Indeed, we observe in the diminution 
of G”, at high frequencies a drop of the line mobility. In 
this case, at constant volume, there is then no alterna¬ 
tive: the curvature of the nano-meniscus must change, 
which we observe in the increase of G’. 

The frequency ranges in our experiments are the ones 
used to operate AFM in dynamical mode. Results pre¬ 
sented here will certainly have to be considered in the 
daily analysis of dynamic mode AFM experiences when 


a water nanobridge is present. Due the importance of 
the observed changes in effective visco-elastic properties 
of water at nanoscale, this conclusion is certainly not lim¬ 
ited to AFM but opens a broader perspective, as capillary 
bridges are now known to be important in surface inter¬ 
actions not only at different scales but also at different 
time scales. 
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